To Jiirg Frohlich whose vision and ideas shaped 
the non-relativistic quantum electrodynamics 



ON RAYLEIGH SCATTERING IN NON-RELATIVISTIC QUANTUM 

ELECTRODYNAMICS 



JEREMY FAUPIN AND ISRAEL MICHAEL SIGAL 

Abstract. In this note we provide details of the proof of the results of |19| (minimal photon velocity bounds 
and - under certain implicit conditions - asymptotic completeness below the ionization threshold, i.e. for 
Rayleigh scattering) for the standard model of non-relativistic quantum electrodynamics, mentioned in that 
paper. 



1. Introduction 

i-C ■ In this note we provide details of the proofs of the main results of [TO] to the standard model of iron- 



relativistic quantum electrodynamics in which particles are minimally coupled to the quantized electromag- 
netic field at energies below the ionization threshold. Recall that in [TO] we proved several lower bounds on 
the growth of the distance of the escaping photons/phonons to the particle system. Using some of these 
results, we proved asymptotic completeness (for Rayleigh scattering) on the states for which the expectation 
of the photon/phonon number is bounded uniformly in time. However, wc provided details only for the 
phonon case. 

Model and assumptions. We consider the standard model of non-relativistic quantum electrodynamics in 
which particles are minimally coupled to the quantized electromagnetic field. The state space for this model 
is given by H := H p ® where H p is the particle state space, F is the bosonic Fock space, T = r(f)) := 
C©^! <8>™f), based on the one-photon space f) := L 2 (R 3 ,C 2 ) ((g)™ stands for the symmetrized tensor product 
of n factors, C 2 accounts for the photon polarization). Its dynamics is generated by the hamiltonian 

n - 

iV aj - KjAtixj)) + U(x) + H f . (1.1) 

j= i i 

Here, rrij and Xj, j = l,...,n, arc the ('bare') particle masses and the particle positions, U(x), x = 
(x\, . . . , x n ), is the total potential affecting the particles, and Kj arc coupling constants related to the 
particle charges. Moreover, := £ * A, where £ is an ultraviolet cut-off satisfying e.g. |(? m £(fc)| < (k)~ 3 , 
\m\ = 0, 1, 2, and A(y) is the quantized vector potential in the Coulomb gauge (div A(y) = 0), describing the 
quantized electromagnetic field and given by 

Mv)--= E / ^=ak)ex{k){e lk -y ax {k) + e- k -yal{k)). 



A=l,2 



Here, u(k) = \k\ denotes the photon dispersion relation (fc is the photon wave vector), A is the polarization, 
and a\(k) and a* x (k) are photon annihilation and creation operators acting on the Fock space T (see Sup- 
plement I for the definition). The operator Hf is the quantum hamiltonian of the quantized electromagnetic 
field, describing the dynamics of the latter, given by Hf := dr(cj), where dr(r) denotes the lifting of a 
one-photon operator r to the photon Fock space, dr(r)|c = for n = and, for n > 1, 

n 

i =1 j—l n—j 

(See Supplement I for definitions and for the expression of dT(r) in terms of the creation and annihilation 
operators. Here and in what follows, the integrals without indication of the domain of integration are taken 
over entire M 3 .) 

l 
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This model goes back to the early days of quantum mechanics (it appears in the review [5T] as a well- 
known model and is elaborated in an important way in [53]); its rigorous analysis was pioneered in |22i 123] 
(see [SU [52] for extensive references) . 

We assume that U(x) € Lf oc (K 3n ) and is either confining or relatively bounded with relative bound w.r.t. 
— A x , so that the particle hamiltonian H p := — Y^Jj—i 2m~^ x j + U{x), an< ^ therefore the total hamiltonian 
H, are self-adjoint. Hence H generates the dynamics through the Schrodinger equation, 

id t iPt = HA- (1-3) 

As initial conditions, "0o, we consider states below the ionization threshold £ (i.e. ipo in the range of the 
spectral projection EV^s) (if)), which is the (largest) energy X > inf cr(H p ) > inf cr(H), below which the 
particle system is well localized: 

||<p)VI*l/(iO||<l, (1.4) 
for any < S < dist(supp /, £) and any / E Cg°((-oo, £)), where (p) := (1 + \p\ 2 ) 1/2 , p := (pi,...,p n ), 
Pj := —iW X j- I n other words, states decay exponentially in the particle coordinates x. To guarantee that 
such an energy exists, we assume that U {x) is such that the particle hamiltonian H p has discrete eigenvalues 
below the essential spectrum ( [35] [6] [7] ) . (£ is close to inf a css (H p ) and is directly defined in [35].) 

We are thus interested in processes like emission and absorption of radiation, or scattering of photons 
on an electron bound by an external potential (created e.g. by an infinitely heavy nucleus), in which the 
particle system (say, an atom or a molecule) is not being ionized. 

Problem. Denote by $^ and Ej the eigenfunctions and the corresponding eigenvalues of the hamiltonian 
H below E, i.e. Ej < E. The following are the key characteristics of the evolution of a physical system, in 
progressive order the refined information they provide and in our context: 

• Local decay stating that some photons are bound to the particle system while others (if any) escape 
to infinity, i.e. the probability that they occupy any bounded region of the physical space tends to 
zero, as t — > oo. 

• Minimal photon velocity bound with speed c stating that, as t — ¥ oo, with probability — > 1, the 
photons are either bound to the particle system or depart from it with the distance > c't, for any 
d < c. 

Similarly, if the probability that at least one photon is at the distance > c"t, c" > c, from the particle 
system vanishes, as t — > oo, we say that the evolution satisfies the maximal photon velocity bound 
with speed c. 

• Asymptotic completeness on the interval (— oo, E) stating that, for any ij; G Ran_E(„ oc x;)(-ff) an d 
any e > 0, there arc photon wave functions fj t <= J-, with a finite number of photons, s.t. the solution 
■ipt = e~ tpo of the Schrodinger equation (| 1 . 3[) satisfies 

limsup||e- itff Vo - $>~^ J % ®- < e- (1-5) 

t— >oo 

3 

In other words, for any e > and with probability > 1 — e, the Schrodinger evolution ip t approaches 
asymptotically a superposition of states in which the particle system with a photon cloud bound to 
it is in one of its bound states, with additional photons (or possibly none) escaping to infinity 
with the velocity of light. 

The reason for e > in (|I.5[) is that, for the state $j <E) S f to be well defined, as one would expect, one 
would have to have a very tight control on the number of photons in /, i.e. the number of photons escaping 
the particle system. (See the remark at the end of Subsection 5.4 of [H] for a more technical explanation.) 
For massive bosons, e > can be dropped (set to zero), as the number of photons can be bound by the 
energy cut-offQ 

We define the photon velocity in terms of its space-time (and sometimes phase-space-time) localization. 
In a quantum theory this is formulated in terms of quantum localization observables and related quantum 
probabilities. We describe the photon position by the operator y := iVfe on L 2 (R 3 ), canonically conjugate 
to the photon momentum k. To test the photon localization, we use the observables dr(ln(y)), where ln(y) 
denotes the characteristic function of a subset il of M 3 . We also use the localization observables T(la(y)), 



For a discussion of scattering of massless bosons in QFT see 
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where T(yJ is the lifting of a one-photon operator \ ( e -g- a smoothed out characteristic function of y) to the 
photon Fock space, defined by 

r(x) = ©~ o(® n x). 

(so that r(e & ) = e dr ^- ) ), and then to the space of the total system. The observablcs dr(ln(y)) and T(la(y)) 
have the following natural properties: 

• dT(l ni un a (l/)) = dT(l ni (l/)) +dT(l na (l/)) and r(l ni (y))r(l 0a (v)) = 0, for fi x and ft 2 disjoint, 

• T g X n (y)T- 1 = X g -i n (y), where T g = r(r 9 ), with t 9 : /(y) -> f{g~ 1 y), and X n (y) stands for either 
dT(l n (y)) or r(l n (y)). 

The observablcs dr(ln(y)) can be interpreted as giving the number of photons in Borel sets f2 C R 3 . They 
are closely related to those used in [351 1331 HZ] (and discussed earlier in [35] and [T]) and are consistent with 
a theoretical description of the detection of photons (usually via the photoelectric effect, see e.g. [50]). The 
quantity T(ln(y))ip) is interpreted as the probability that the photons are in the set Q in the state vj). 
This said, we should mention that the subject of photon localization is still far from being settled]! 

The fact that for photons the observables we use depend on the choice of polarization vector fields, £\(k), 
A = 1, 20 is not an impediment here as our results imply analogous results for e.g. localization observables of 
Mandel [49] and of Amrein and Jauch and Piron [TJEl]: dr(/^ an ) and dr(/^ jp ), where f$ an := PHq^P 1 - 
and /^ JP := ln(y) n P^, respectively, acting in the Fock space based on the space () = L 2 ransv (R 3 ; C 3 ) := 
{/ e L 2 (R 3 ;C 3 ) : k ■ f(k) = 0} instead of f) = L 2 (R 3 ;C 2 ). Here P 1 ^ : f(k) -> f(k) - \k\~ 2 k k ■ f(k) is 
the orthogonal projection on the transverse vector fields and, for two orthogonal projections P± and P2, the 
symbol P\ n P2 stand for the orthogonal projection on the largest subspace contained in Ran P\ and Ran P2 . 

We say that the system obeys the minimal photon velocity bound if the Schrddingcr evolution, ip t = 
e~ lt -ff-0O! beys the estimates 

/oo 
dt t- Q '\\dr( x ^ =1 )h t \\ 2 <\\u\l 

for some norm HV'olloi some < a' < 1, and for any a > and c > such that either a < 1, or a = 1 and 
c < 1. In other words there are no photons which either diffuse or propagate with speed < 1. Here Xfi( v ) 
denotes a smoothed out characteristic function of the set f2, which is defined at the end of the introduction. 
The maximal velocity estimate, as proven in [10j . states that, for any d > 1, and 7 < min(g(l — 

||dr( X M >1 )^ f || < rT||(dr« y » + 1)^ ||. 

c' t — 

Considerable progress has been made in understanding the asymptotic dynamics of non-rclativistic particle 
systems coupled to quantized electromagnetic or phonon field. The local decay property was proven in 
[HI [HI HH UHl [3D [32> by the combination of the renormalization group and positive commutator 
methods. The maximal velocity estimate was proven in |10j . 

An important breakthrough was achieved recently in |14j . where the authors proved relaxation to the 
ground state and uniform bounds on the number of emitted massless bosons in the spin-boson model. 
(Importance of both questions was emphasized earlier by Jiirg Frdhlich.) 

In quantum field theory, asymptotic completeness was proven for (a small perturbation of) a solvable 
model involving a harmonic oscillator (see [31 EI]), and for models involving massive boson fields ( jTTJ [2S1 
[26l[27]). Moreover, [33] obtained some important results for massless bosons (the Nelson model) in confined 
potentials (see below for a more detailed discussion). Motivated by the many-body quantum scattering, 
[T71 [3SJ [23 [13 [33] defined the main notions of scattering theory on Fock spaces, such as wave operators, 
asymptotic completeness and propagation estimates. 

Results. Now we formulate our results. It is known (see [313S]) that the operator H has a unique ground 
state (denoted here as 'I'gs) and that generically (e.g. under the Fermi Golden Rule condition) it has no 



The issue of localizability of photons is a tricky one and has been intensely discussed in the literature since the 1930 and 
1932 papers by Landau and Peierls [46] and Pauli |52| (see also a review in |45j1. A set of axioms for localization observablcs 
was proposed by Newton and Wigner 1511 and Wightman 1641 and further generalized by Jauch and Piron 1441 . Observables 
describing localization of massless particles, satisfying the Jauch-Piron version of the Wightman axioms, were constructed by 
Amrein in [l]. 

3 Since polarization vector fields are not smooth, using them to reduce the results from one set of localization observables to 
another would limit the possible time decay. However, these vector fields can be avoided by using the approach of |48j . 
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eigenvalues in the interval (E gs ,T), where E gs is the ground state energy (see [SJ [25J [33] ) . We assume that 
this is exactly the case: 

Fermi's Golden Rule ([3(7]) holds. (1.6) 
(If the particle system has an infinite number of eigenvalues accumulating to its ionization threshold - which 
is the bottom of its essential spectrum - then, to rule out the eigenvalues in the spectral interval of interest, 
we should replace S by S — e for some fixed e. This is understood from now on.) 
We will consider the following sets of initial conditions 

Ti := {Vo G /(^^(dr^- 1 )*), for some / e Cg°((-oo,E))}, 

and 

T 2 := {^o € f(H)D(dT((y))), for some / £ C%>((E gs , Ti))}. 

For A > —C, let HV'olU := \\(A + C + 1)^i/}q\\. We define v{p) > as the smallest real number satisfying 
the inequality 

^ u dV{Lo^ t )<t^\\M\% (1.7) 
for any ipo £ Rani5(_ oo S ) (H), where \\ip\\ P := HV'lldr(wP)- It was shown in pjj] (see also (A.l) of Appendix 
A of [19]) that 

«A < ^ 

(this generalizes an earlier bound due to [33]). Also, the bound ||V't||i?, ^ HVtoHff shows that (|1.7p holds for 
p = 1 with i/(l) = 0. With v(p) defined by (jTTTJ) , we prove the following two results. 

Theorem 1.1 (Minimal photon velocity bound). Let either j3 = 1 and c < 1, or 

i(5 + K-l)-K0)) </?<!• (1.8) 
6 

Then for any initial condition ipo £ T±, the Schrddinger evolution, ip t , satisfies, for any a > 1, the following 
estimate 

f>OC 

dt t-?- av f)||dT(x i*i - x )Htf < HoW-i- (1.9) 

Theorem 1.2 (Weak minimal photon escape velocity estimate). Let the coupling constants Kj be sufficiently 
small. Assume (|1.6p and u(—l) < a < 1 — i>(0). Then for any initial condition tf>Q £ T2, the Schrddinger 
evolution, tp t , satisfies the estimate 

llrtx^JVtH < t-< {\\Mlr {{ y)) + II^HSr (5)a ), (1.10) 

where 7 < | min(l — a — u(0), |(a — z/(0) — v (— 1))) and b := ■ y + y ■ k). 
Remarks. 

1) The estimates ([T9]l and (fTTTOl) are sharp if i/(0) = and ^(-1) = 0. 

2) The estimate (jl.101) states that, as t — s- 00, with probability — > 1, either all photons are attached to 
the particle system in the ground state, or at least one photon departs the particle system with the distance 
growing at least as 0{t a ). 

3) The norm on the r.h.s. of (|1.10j) above is different (weaker) than the one given in [TH], however 
the proof can be easily upgraded to the weaker norm indicated, by omitting the unnecessary estimate 
IIV>o||^ )2 < IIV-oll^c^ for Ifo e D{H). 

Let N := dr(l) be the photon number operator and K := dr(a> -1 ) be the photon low momentum number 
operator. Our next result is 



Theorem 1.3 (Asymptotic Completeness). Let the coupling constants Kj be sufficiently small. A 



ssume 



(|1.6|) and suppose that, uniformly in t £ [0,oo), either 

\\n 1/2 M<\\n 1/2 M + UoI (1.11) 

for any Vo G f^DiN 1 / 2 ), with f £ C^((E gs , 53)), or 

I|A-^|I<1, (1.12) 
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for any ipo G T>, where T> is such that DO D(dr(u 1 / 2 (y)uj 1 ^ 2 )^ 2 ) is dense in Ran i/r-oo Then the 
asymptotic completeness holds on RanS(_ 00] s)(iJ). 

The advantage of Assumption (|1.12p is that the uniform bound on K = dr(w _1 ) is required to hold 
only for an arbitrary dense set of initial states and, as a result, can be verified for the massless spin-boson 
model by modifying slightly the proof of [14] (see [20]). Hence the asymptotic completeness for the massless 
spin-boson model holds with no implicit conditions. 

As we see from the results above, the uniform bounds, or (|1.12|) . on the number of photons emerge as 

the remaining stumbling blocks to proving the asymptotic completeness without qualifications. The difficulty 
in proving these bounds for massless fields is due to the same infrared problem which pervades this field and 
which was successfully tackled in other central issues, such as the theory of ground states and resonances 
(see [5] [56] for reviews), the local decay and the maximal velocity bound. 

For massive bosons, the inequality (jl.llj) (as well as (|1.7j) . with z/(0) = 0) is easily proven and the proof 
below simplifies considerably as well. In this case, the result is unconditional. It was first proven in |17j for 
models with confined particles, and in |25) for Rayleigh scattering. 

Comparison with earlier results. For models involving massive bosons fields, some minimal velocity 
estimates are proven in |17j . For massless bosons, Theorems 11.11 and 11.21 seem to be new. As was mentioned 
above, asymptotic completeness was proven for (a small perturbation of) a solvable model involving a 
harmonic oscillator (see [31 I61j). and, for models involving massive boson fields, in |17] for confined systems, 
in |25] below the ionization threshold for non-confined systems, and in |26j for Compton scattering. 

The paper [33J treats the Nelson model given by the state space H := "K v ® J 7 , where Hp := L 2 (R 3 ), T is 
the bosonic Fock space based on the one-phonon space f) := L 2 (M 3 ,C), and the hamiltonian 

H := H p + Hf + 1(g), (1.13) 

acting on H. Here H p := ~-^A x + V(x), with V(x) growing at infinity as V(x) > Co|x| 2a — C\, Cq > 0, 
a > 0, acting on L 2 (R 3 ), Hf = dT(uj), where uj(k) = \k\ and 1(g) is given by 

1(9) ■= J(g*(k)®a(k)+g(k)®a*(k))dk, (1.14) 

where a*(k) and a(k) are the phonon creation and annihilation operators acting on J 7 , with abstract condi- 
tions on the coupling function g (allowing a coupling function of the form g(k) = |fc|^£(fc)e* fc:i; , where £(fc) 
is the ultraviolet cut-off, with various conditions on p depending on the results involved). In this case, in 
particular, the ionization threshold £ is equal to oo. 

We reproduce the main results of [33] (Theorems 12.4, 12.5 and 13.3), which are coached in different 
terms than ours and present another important view of the subject. Let /, /o £ C°°(K) such that < 
f,fo < 1, /' > 0, / = for a < a , f = 1 for a > a u f^< 0, f = 1 for s < a u f = for 
s > a 2 , with < a < a\ < a 2 . Let P+ := inf c<c 'P+, with P+ := s-lim e ^ e' 1 R+(e~ 1 ), Rti^ 1 ) ■= 
s-Hm^oo e ttH (B ct + \)~^e~ ltH , B ct := dr(6 ct ), b ct := /(^) and r+(/ ) := s-lim^ e UH T(f , c ,, t ) e - ltH , 
where /o, c ',t '■— fo( 7p° * )■ Then Proposition 12.2 and Theorem 12.3 state that the operators P c + exist 
provided p > — ^j, are independent of the choice of /, and are orthogonal projections commuting with H. 
Furthermore, let JC + := {$ € W,a±(h)& = 0, V7i € i)} (called in [33] the set of asymptotic vacua), where 
(formally) a±(h) := s-\va\ t ->±oo e. ltH a(e" ltul h)e~ ltH and := RanP+ (the spaces containing states with 
only a finite number of photons in the region {\y\ > c't} as t — > oo, for all d > c). Assuming a > 1 and /i > 0, 
Theorems 12.4 and 12.5 state that the operator Tj(/o) exists and is equal to the orthogonal projection on 
the space /C+ := JC + n provided < c < c' < 1 and p > j^ff- Assuming in addition that the Mourre 
estimate 1a(H)[H, iB]l&(H) > cq1a(H) + R holds on an open interval Act, with the conjugate operator 
B := dr(6), Co > and R £ % compact, then for < c < c(A,co), one has lj^(H)K^ = 1a(^) - H pp , 
where "H pp is the pure point spectrum eigenspace of H. (The latter property is called in |33j geometric 
asymptotic completeness. Combining results of [8] [28] one can probably prove a Mourre estimate, with B 
as conjugate operator, in any spectral interval above E gs and below £ and for the coupling function g given 
by g(k) = \k\^£(k)e lkx , with p > 1/2.) 

Our approach is similar to the one of |33j in as much as it also originates in ideas of the quantum 
many-body scattering theory. At this the similarities end. 
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Approach and organization of the paper. In [TI5] we gave the details of the proofs of similar results for 
the hamiltonian (|1.13[) - ([1.14[) . with the coupling operators g(k) satisfying 

\\v H d a g(k)\\n P < \kr^(k)- 2 -^ \a\ < 2, (1.15) 

where /i > —1/2 or /i > 0, and rj is an estimating operator on the particle space "Hp (a bounded, positive 
operator with unbounded inverse), satisfying 

\\r n f{H)\\ <i, 

for any n = 1, 2 and / S C§°((— oo, £)). As was mentioned in [19], to extend these proofs to the hamiltonian 
we use, similarly to [TU], the generalized Pauli-Fierz canonical transform to map the hamiltonian fll.l[) 
to a hamiltonian H with milder infrared behaviour (Section [2]). Then we extend readily techniques of |19j to 
H and relate the results for H with those for (Section [3]). We also include Supplement I defining and 

discussing the creation and annihilation operators and presenting some standard commutation relations. 

Notations. For functions A and B, we will use the notation A < B signifying that A < CB for some 
absolute (numerical) constant < C < oo. The symbol Ea stands for the characteristic function of a set 
A, while x <i denotes a smoothed out characteristic function of the interval (— oo, 1], that is it is in C°°(K), 
is non-decreasing, and = 1 if x < 1/2 and = if x > 1. Moreover, x >i : = 1 — X<i an d X =i stands for 
the derivative of x >i- Given a self-adjoint operator a and a real number a, we write Xa<a := X-<ii an d 
likewise for Xa>a- Finally, D(A) denotes the domain of an operator A. 



2. Generalized Pauli-Fierz transformation 

We consider for simplicity a single negatively charged particle in an external potential. Then, absorbing 
the absolute value of the particle charge into the vector potential A(x) := A^(x) and choosing units such 
that the electron mass is m = 1/2, the hamiltonian (jl.lj) becomes 

H={p + A(x)) 2 + H f + U{x). (2.1) 

The coupling function g^ ed (k,X) := \k\^ 1 ^ 2 £ i (k)e\(k)e lk ' x in (|2.1I) is more singular in the infrared than is 
allowed by our techniques (/i > 0). To go around this problem we use the (unitary) generalized Pauli-Fierz 
transformation (see |55j ) 

U:=e^ {q *\ (2.2) 

to pass from H, given in (|2.ip . to the new hamiltonian H := UHU*, where $(/i) is the operator- valued 
field, $(/i) := -^(a*(h) + a(h)) and the function q x (k,X) is defined as follows. Let tp € C°°(]R;IR) be a 
non-decreasing function such that ip(r) = r if \r\ < 1/2 and \<p(r)\ = 1 if |r| > 1. For < u < 1/2, we define 

q x ( k ,\):=-^-<p(\kre x (k)-x). (2.3) 

|fc| 2 +,y 

We note that the definition of <f>(/i) gives A(x) = $(g§ od ). Using p.l[) and p.2[) . we compute 

H=(p + A{x)f + E(x) + H f + V(x), (2.4) 

where 

A{x) := A) := 5 2 cd (fc, A) - V^^fe, A), 

E(x) := $(e K ), e x (A;, A) := i|fc|g K (A:, A), 

:= U{x) + | Ea=i,2 Jr3 A)| 2 dfc. 

The operator if is self-adjoint with domain V(H) = V(H) = V{p 2 + H f ) (see [39l [40] V 

Now, the coupling functions (form factors) g x {k, A) and e x (k, A) in the transformed hamiltonian, H, satisfy 
the estimates that are better behaved in the infrared (|10|): 

\drgx(k,\)\<(k)- 3 \k\?-W(x)l+\ m \, 

\dZ l e x (k,\)\<(kr 3 \k\?-W(x) 1+ W. 
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We see that the new hamiltonian ([1I4"]) is of the form (|1.13l) . H = H p + H f + 1(g), with H p := -A + V(x), 
H f = dr(w) and 

Kg) -=J2 d k(i)dk' u) g l:j (k il} ,!i {j} ) ® a*(fc (4) )a(fc' b) ), 

where the summation in i, j ranges over the set > 0, 1 < i + J < 2, fc^ p j := . . . , fc p ), fc^ := (kj, Xj), 
I d k {p ) '■= WiY^xJ dk ii a *(k( p) ) '■= Y\ P i a *{kj) if P > land= 1, if p = 0, a#{k } ) := a*(k 3 ), and g := {g tJ ). 
The coupling operators, = gij{k(i)-,k(j)) obey 

9ij (k(i) > k(j) ) = 9ji (k{ 3 ) i k{i) ) ! 

and satisfy the estimates 

\\{p) i+j - 2 {x)-^d a gij (k ii+j) )\\n P < E n(NW- 3 )iM" w , (2.5) 

m=l 1=1 

with /i = 1/2, \a\ < 2, and 1 < i + j < 2. 

The bound (TO]) holds for both (f2TTj) and (|2~4"|) . With (fl~4]) and (|2~5|) . it is easy to extend Theorems O 
and 11.21 to the operator H and then translate the result to (|2.1[) . Similarly one proceeds with Theorem ll.3l 
We explain all this in Section [3] 



3. Proof of Theorems ll.lhfl~3l 

3.1. Proof of Theorems [TTTlfPl Using ([23]) . it is not difficult to extend Theorems O and O to the 
operator H (cf. [TU], where one can find many needed, or similar, estimates). Then we translate Theorems 
l.ll and ri.2l from H to the QED hamilonian (|2.ip . by using the following estimates ([10]) 



dr( Xl (v))H' 

T( X 2(v))h 



< 



+ t- ad \m 2 , 



< 



u^,v{ X2 (v))U4>) + r ad u\\ 2 



(3.1) 
(3.2) 



where v 



valid for any functions Xii v ) an d X2( v ) supported in {\v\ < e} and {\v\ > e}, respectively, 



for some e > 0, for any i/j G f(H)U, with / 6 C§°((-oo, £)), and for < d < 1/2. (J3HJ follows from 
estimates of Section 2 of [TU] and (|3.2[) can be obtained similarly (see (|I.2[) and (|I.3[) ). Using these estimates 
for ip t = e~ i/)q, with an initial condition ipo in either Ti or T2, together with Ue~ i/jq = e~ Utpo, and 
applying Theorems 11.11 and 11.21 for H to the first terms on the r.h.s., we see that, to obtain Theorems 11.11 
and 11.21 for the hamiltonian (|2.1[) . we need, in addition, the estimates 



^,W*dr(o;- 1 )WV) < (if>, (dr(o;- 1 
4',U*dT((y))U^) < (^(dT«y))- 
W*dT(6)WV|| £ || (dT(6) + (a;) 2 )^! 



1)V>, 
c> 2 )^), 



(3.3) 
(3.4) 
(3.5) 



where, recall, & = ^(k ■ y + y ■ k). 



To prove p.3[) , we see that, by (|I.2[) . we have 

W*dr(o;- 1 )W = e^^-Jdrfw-^e - **^ 5 = dr^ 1 ) - $(^Vt) 



I^VII 2 , 



(Since w -1 ^ ^ f), the field operator ^(iuj~ 1 q x ) is not well-defined and therefore this formula should be 
modified by introducing, for instance, an infrared cutoff parameter a into q x . One then removes it at the 
end of the estimates. Since such a procedure is standard, we omit it here.) This relation, together with 



< 



,-3-2i/+s/ 



(k)- e dk -||di>^)^||IMI, 



for any e > 0, which follows from the bounds of Lemma 1.1 of Supplement I, and 



\u 2 fell ft 



< 



'(A)- 



(3.6) 
(3.7) 
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implies p.3[) . 

To prove (|3.4[) and p.5[) , we proceed similarly, using, instead of (|3.6[) and (|3.7[) . the estimates 



K^*(i(»)&)V)| < ( / w- 2 - 2 "<fc)- 8 dfc)*||dr(w- 1 )m|||<x)^|| 

< ( f uj- 2 - 2 »(k)- 6 dkf\\dT((y))h\\\\(xW\\, 



and 

H*(&feMI ;S ( / ^ 2 -^(fc)" 6 *)'ll(^>(^/ + i)^ll, 

3.2. Proof of Theorem 11.31 We present the parts of the proof of Theorem 11.31 for the hamiltonian (|2.1[) 
which differ from that for the hamiltonian (|1.13[ l- (|1.14[) . To begin with, the existence and the properties of 
the Deift-Simon wave operators, 

W± := s-lim W(i), with W(t) := e u "t (j)e" 4iff , 

t— >±oo 

on Ran^oo where H = H (g) 1 + 1(g) 77/ and the operators f and j = (jo,joo) are defined in Subsection 

5.1 of [19], are equivalent to the existence and the properties of the modified Dcift-Simon wave operators 

4 mod) := s-lim (e"**^ ® l) e ^f (j>- 4 "V*^\ (3.8) 
on Ran { _ OCiS) ( J ff) (where H = er^^He 1 ^^ is given in $£ty). 



To prove the existence of W± mod ' 1 , we observe that, due to 

f (j)a*(h) = (a*(j h) 1 + 1® o^O'ooh))^*). ( 3 - 9 ) 
where a# stands for a or a*, we have T(j)&(h) = ^(h)T(j), where 

Hh) := <P{j h) <g> 1 + 1 ® $(ioo/i), 

which in turn implies that 

f (j) e i$ W = e^Wf (j). 

Therefore 



e ^ (mod) f(j) e -^ +Rem t , (3.10) 



where f/( mod ) := If <g) 1 + 1 g) iJf and Rcm t := (er^^ <g> l)e itA (e i4, (^) - e i$ ( fo ) ® l)f {j)e- ltil . We claim 
that 

s-lim Rcmt = 0. (3-11) 

t— >±oo 

Indeed, let R := $(g x ) - <g> 1 = $((j - ® 1 + 1 ® ^(joofe) and TV := TV ® 1 + 1 ® N. Using ([13]), 

Lemma 1.1 and Lemma B.6 of [115] (see also [TU] Lemma 3.1]), we obtain 

\\R(N + 1)-*|| < HOo - + lliocfellf, < t- QT (*) 1+T , 

for any r < 1. From this estimate and the relation e i4 ^) - e**^ 1 = -i /J dse^ -8 ^^ R(e si *^ ® 1), 
it is not difficult to deduce that 

||( e **(9-) - e I$fe) ® l)(/V+ (x) 2+2r + I)" 1 )) < i~ Qr . 

Furthermore, we have (TV + (a;) 2 + 2r + l)f (j) = f(j)(N + (x) 2+2t + 1), and, as in Corollary A.3 of [15], 
one can verify that ||iVe- lt ^V'o|| < t 2/5 |K'o||-i for any Vo e /(# )-D(# 1/2 ), / e Cg°((-oo,S)). Using in 
addition that |j (x) 2+2r /(i?)|| < oo, it follows that Rem t strongly converges to on Rffli(_ m j;)(-ff) provided 
that err > 2/5. 
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The equations ^E\i, (|3~TU]) and (|3~TT]) imply 



(mod) 



s-lim e 

t— J-±oo 



; t ^(mod) 



'± T(j)e- UM . (3.12) 

The proofs of the existence and properties of the Dcift-Simon wave operators (|3.12p are then done as for W± 
in [19] (see Theorem 5.1 of [H]), but with Gi, defined in (5.16) of [19], replaced by 

Gi := (1(g) ® l)f(j) - t(j)I(g), (3.13) 

where 1(g) := p ■ A(x) + A(x) ■ p + A(x) 2 + E(x) is coming from the decomposition 

H = H p + H f + I(g). 

to derive 

f (j)I(g) = i(g)f( 3 ), 



We use 



where 

1(g) := p ■ A(x) + A(x) ■ p + A(xf + E(x), 

with A(x) := ®(jog x ) <g) 1 + 1 <g) ^OocEfa), E(x) := $0' e K ) ® 1 + 1 ® ^C^e^). This together with (|3~T5]l 
gives 

G 1 =7 # ( fl )f(i) ! 

where 

:= P ■ it#(a:) + • p + A # (x) ■ A(x) + A(x) ■ A # (x) + E # (x), 

with A # (x) := $((1 - j )g x ) ® 1 - 1 ® ^(joo9x), E # (x) := $((1 - joK) <g) 1 - 1 (8 ^(jooe*). Using the last 
expression, Lemma 1.1 and the estimates 

\\joog x h<t- ar (x)i+^ \\(i-jo)g x h<t- aT (x)i + \ 

Hiooexllft <t- ar (x) 1+ \ \\(l-Jo)e x h<t~ aT (x) 1+ \ 
with r < 2 (see Lemma B.6 of [19]), we arrive at the following refinement of estimate (5.31) of [19] 

\\f(H)G 1 (N + l)- 1/2 \\ <t- aT . 

The remainder of the proof goes the same way as the proof of Theorem 5.1 of [19]. Hence W^ mod ' ) and 
therefore W± exist. 

Now, we comment on the proof of the key Theorem 5.4 of [19] in the QED case. It goes in the same way as 
for the phonon model (|1.13[) - (|1.15[) (see [IH]), until the point where we have to show that ||dL(6^) 1 / 2 P gs || = 
0(t K ) in the present case. This estimate can be proven by using the generalized Pauli-Fierz transformation 
together with (|I.3[) . to obtain 

1 



|dT(&?)** 



$ gs , (dT(b 2 ) - ®(ib 2 e q x ) + ^fe, 



(3.14) 



where $ gs := UQ gs . Using Lemma 1.1 of Supplement I, (|1.4[) and the fact that $ gs £ D(N 1 ^ 2 ), we can 
estimate the second term of the r.h.s. of p,14j) as 

($ gs ,$(^<7 x )$ gs }| < ||( a; ) 3 $ ff ||||( a: )- a $(i^ fe )(JV + l)-4||||(JV- + l)3$ g|1 || <t 2K . 
Likewise, |($ gs , (b 2 q x ,q x )tj^ gs }\ < t 2ft . To estimate the first term of the r.h.s. of (|3.14[) , we write 



\dT(b 2 e )^ ss \\ 2 =Y. I ||MA(fc)$ gs || 2 rffc. 



Applying the standard pull-through formula gives 

a A (fc)$ gs = -(H-E gs + \k\y 1 ((p + A(x))-g x (k) + e x (k))^ gs . 
We then easily conclude that ||& e aA(fc)$ gs || f, = 0(t K ) in the same way as in Corollary C.4 of 
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Supplement I. Photon creation and annihilation operators 

Recall that the propagation speed of the light and the Planck constant divided by 27r are set equal to 1. 
We use the momentum representation and write functions from the space f) = L 2 (M. 3 , C 2 ) as u(k, A), where 
k £ M. 3 is the wave vector or momentum of the photon and A £ { — 1, +1} is its polarization. 

With each function / el), one associates creation and annihilation operators a(f) and &*(/) defined, 
for u £ ®"f), as 

d*(/):ii->\m + l/®s« and a(f) : u — > \/n(f,u)t ! , 
with (/, u)tf := ^2\ =1 2 J dk f(k,X)u n (k,X,kx, X±, . . . ,k n —i, X n —x). They are unbounded, densely defined 
operators of r(h), adjoint of each other (with respect to the natural scalar product in JF) and satisfy the 
canonical commutation relations (CCR): 

[a*(f),a*(g)] = 0, [a(f),a*(g)} = (f,g), 
where d*=aora*. Since a(/) is anti-linear and a*(f) is linear in /, we write formally 

«(/)=E / dkf^X)a x (k), a*(f)=Y^ fdkf(k,X)a* x (k), 

\=l,2 \=1,2 

for photons. Here a x (k) and a* x (k) are unbounded, operator- valued distributions, which obey (again formally) 
the canonical commutation relations: 

[a*(k),a*(k')] =0, [a x (k),a* x ,(k')] = 5 x ,y5(k - k'), 

where = a x or a* x . 

Given an operator r acting on the one-particle space f), the operator dr(r) defined on the Fock space 
T by (|1.2[) can be written (formally) as dT(r) := ^2 X=12 J dk a* x (k)Ta x (k). Here the operator r acts on 
the (k, A)-variable (dr(r) is the second quantization of r). The precise meaning of the latter expression is 
(|1.2[) . In particular, one can rewrite the quantum Hamiltonian Hf in terms of the creation and annihilation 
operators, a and a*, as 

H f =^2 I dka* x (k)uj(k)ax(k). 

A=l,2*' 

Commutators of two dr operators reduces to commutators of the one-particle operators: 

[dr(r),dr(r')] =dr([r,r']). 

Let r be a one-photon self-adjoint operator. The following commutation relations involving the field 
operator $(/) = Tfe(o*(f) + a (/)) can be readily derived from the definitions of the operators involved: 

[$(/), <&(.<?)] =ilm(f,9)i, 
[$(/), dr(r)] = **(«-/), 

[r(r), $(/)] = r(r)a((l - r)f) a*((l - r)f)T(r). 
Exponentiating these relations, we obtain 

e^<S>(g)e-^ = <%) - Im(f, g) h , (LI) 

e l *(«dr(r)e- J *( / ) = dT(r) - §(ir/) + \ Re(r/, /)„ (1.2) 

e »*(/)r(r)e- i *( /) - F(r) + f dse ls<I,(/) (r(r)a((l - r)f) - o*((l - T)f)T{T))e~ sl ^ . (1.3) 

Jo 

Finally, we have the following standard estimates for annihilation and creation operators a(/) and a*(f), 
whose proof can be found, for instance, in [7] , [351 Section 3] , [35] : 

Lemma 1.1 For any f £ f) such that uj~ p / 2 f £ f), the operators a*(/)(dr(o; p ) + l) -1 / 2 , where a#(/) stands 
for a*(f) or a(f), extend to bounded operators on % satisfying 

||a(/)(dlV) + l)-*|| < Hw-^/llft. 

|| a * (/)(drK) + i)-i II <|| w -"/ 2 /lk + 11/11.,. 
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//, in addition, g € f) is such that uj p / 2 g 6 t), the operators a# (f)a& (g)(dT(uj p ) + 1) 1 extend to bounded 
operators on % satisfying 

llaC/MsXdlVO + l)- 1 !! < Wu-^fUu-^gU, 

\\a*(f)a(g)(dT(u>») + 1)^ < (ll^Vlk + Wfh) W^^gh, 

\\a*(f)a*(g)(dT(coP) + l)^ < (||^/ 2 /|| t, + Wfh) (W^Qh + Uh)- 
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